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Abstract 

The  stress  field  in  the  vicinity  of  a  broken  fiber-reinforced  composite  is 
analyzed  by  means  of  a  shear-lag  model.  The  broken  filament  is  posi¬ 
tioned  eccentrically  relative  to  its  neighboring  fibers  to  simulate  the  com¬ 
monplace  non-uniformity  of  fiber  spacing  within  the  transverse  plane.  It 
is  shown  that  a  fiber  break  gives  rise  to  severe  bending,  in  addition  to 
tension,  in  the  neighboring  fibers  -  with  a  substantial  overstress  focused 
on  the  nearest  unbroken  filament.  The  complex  nature  of  the  stress  field, 
which  is  caused  by  the  failure  of  a  fiber  within  a  composite  casts  doubt  on 
the  applicability  of  failure  statistics  derived  from  tensile  failure  data  of 
single  fibers,  which  is  commonly  used  to  predict  the  strength  of  compos¬ 
ites. 


*Pcmiaiicnt  Address:  Holon  Institute  of  Technological  Education.  Affiliated  with  Tcl- 
Aviv  University,  P.  O.  Box  303,  Holon  38102,  ISRAEL. 


INI'RQDUCTION 


Concentric  shear-lag  models  have  been  used  extensively  in  the  analysis  of 
composites,  in  particular,  in  the  context  of  micromechanical  aspects  of  fail¬ 
ure  (see,  for  example,  Hull,  1981).  Indeed,  the  interplay  between  mi¬ 
cromechanics  and  the  statistics  of  fiber  failure  appears  to  be  the  most  dif¬ 
ficult  aspect  in  modeling  the  strength  of  fibrous  composites. 

Due  to  the  complexity  of  the  phenomenon,  the  theoretical  analysis  invokes 
a  variety  of  assumptions,  the  influence  of  which  may  be  difficult  to  assess. 
Oiiv  of  tlie  basic  appioachcs  was  set  forth  by  Giicer  and  Gurland  (1962) 
who  represented  a  material  as  a  series  of  layers.  Each  layer  was  consid¬ 
ered  to  consist  of  fiber-like  elements  loaded  in  parallel,  which  allows  for 
the  use  of  the  bundle  theory  (Daniels,  1945)  to  analyze  its  strength.  The 
overall  failure  of  the  layered  structure  may  then  be  predicted  by  the 
weakest  link  approach.  This  theory  has  left,  however,  the  layer  thickness 
unspecified.  Rosen  (1965),  proposed  to  identify  this  parameter  as  the  so- 
called  fiber  ineffective  length,  8,  and  provided  a  definition  of  this  quantity. 

Since  then  this  approach  which  combined  micromechanical  and  statistical 
considerations,  has  been  of  particular  interest  in  failure  investigations. 
Garg  et,  al  (1973),  provided  a  good  exposition  of  the  approach,  which  re¬ 
mained  of  particular  interest  up  to  the  present  time,  such  as  in  the  recent 
investigations  of  composite  failure  by  Prewo  (1986;  and  Schweifert  and 
Steif  (1990). 

As  noted  earlier,  the  models  of  statistical  strength  are  essentially  approxi¬ 
mate  and  invoke  various  assumptions.  A  fortuitous  agreement  with  exper¬ 
iments  may  therefore  well  occur.  It  is  instructive  that  a  prediction  of  fail¬ 
ure  which  is  based  upon  the  simple-minded  rule  of  mixtures  appears  to  fit 
the  experiments  reported  by  Prewo  (1986)  fairly  well,  while  a  much  more 
sophisticated  approach  by  Zweben  and  Rosen  (1970)  provides  good  results 
only  if  the  ineffective  length  is  taken  to  be  orders  of  magnitude  larger  than 
the  value  derived  from  the  theoretical  analysis.  Thus,  it  becomes  clear 
that  there  may  be  important  micromechanical  and  other  factors,  which 
were  overlooked  by  the  analysis. 

To  discuss  these  factors  we  recall  the  basic  result  given  by  Rosen  (1965)  as 
follows 

^  '/o  /  r  (u  )  (1) 
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Here  o*  and  o  are  the  statistical  mode  of  composite  failure  stress  and 
the  mean  tensile  strength  of  an  individual  fiber,  respectively,  S  and  L  are 

the  ineffective  length  and  the  fiber  length,  respectively,  P  is  an  empirical 

Weibull  parameter,  and  r(  ■ )  is  the  Gamma  function.  One  observes  from 
Equation  (1)  that  the  composite  mechanical  behavior  manifests  itself 
through  the  ineffective  length,  S,  while  the  statistical  nature  of  the  phe¬ 
nomenon  is  introduced  through  the  Weibull  parameter,  which  is  ob¬ 
tained  from  failure  data  for  "dry"  fibers  in  pure  tension. 

There  appear  to  be  two  major  arguments  which  cast  doubt  on  the  accuracy 
of  Equation  (1)  and  its  underlying  methodology.  The  first  argument,  for¬ 
warded  by  McCarthy  and  Orringer  (1975),  indicates  that  statistical 
strength  data  obtained  for  tensile  loading  of  fibers  involve  samples  with 
lengths  of  about  one  inch.  On  the  other  hand,  according  to  the  analysis  of 
Rosen  (1965),  the  ineffective  length,  8,  which  specifies  the  layer  thickness 
in  Equation  (1),  is  usually  between  10  and  40  microns.  This  brings  into 

question  the  validity  of  the  value  of  P ,  as  obtained  from  fiber  tensile  tests, 
in  predicting  the  strength  of  composites. 

The  second  argument  deals  with  the  absence  of  load  concentration  effects 
in  Equation  (1).  Zweben  (1968)  and  Zweben  and  Rosen  (1970),  attempted 
to  modify  the  theory  in  order  to  account  for  these  effects.  While  the  over¬ 
load  suffered  by  the  fibers  adjacent  to  a  broken  filament  does  indeed  play 
an  essential  role  in  the  failure  mechanism,  the  above  works  still  do  not 
provide  a  satisfactory  agreement  with  experiments  (Prewo,  1986). 

The  above  observations  show  that  there  may  be  factors  which  are  over¬ 
looked  by  the  previous  analyses.  Among  these  we  would  particularly  note 
the  randomness  of  the  fiber  location  in  the  transverse  plane.  Typical  mi¬ 
crographs,  for  graphite/BMI  and  graphite/epoxy  composites  are  shown  in 
Figures  (la)  and  (lb),  while  results  of  Monte-Carlo  simulations  are  exhib¬ 
ited  in  Figures  (2a)  and  (2b).  It  is  obvious  that  since  fibers  are  not  uni¬ 
formly  spaced  within  the  cross  sections,  the  concentric  shear-lag  model 
does  not  account  for  this  commonplace  phenomenon  within  composites. 

The  eccentric  model,  introduced  herein,  reveals  essential  effects  induced 
by  the  above  mentioned  randomness.  The  present  investigation  indicates 
that  in  contrast  with  the  premises  of  bundle  theory,  the  failure  of  a  fila¬ 
ment  in  a  composite  causes  a  complicated  stress  field  in  the  vicinity  of 
break.  Specifically,  the  neighboring  fibers  suffer  extensive  bending,  in 
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addition  to  tension,  and  may  therefore  be  subjected  to  failure  mechanisms 
other  than  thc'^o  associated  with  pure  tension.  Though  bending  in  neigh¬ 
boring  fibers  occurs  even  in  concentric  configurations  around  a  broken  fil¬ 
ament,  eccentricity  serves  to  accentuate  bending  effects.  It  is  shown  in  the 
present  paper  that,  with  rather  commonplace  eccentricities,  the  overload 
due  to  bending  may  be  quite  substantial,  with  values  of  about  twice  those 
of  the  intact  case. 

Consequently,  the  present  micromechanical  analysis  casts  further  doubt  on 
the  applicability  of  the  commonly  used  statistical  data  obtained  for  free 
fibers  in  pure  tension  to  the  prediction  of  strength  of  fibrous  composites 

2.  ANALYSIS 

Shear-lag  models,  which  abound  in  the  mechanics  literature  provide  useful 
approximations  for  the  mechanism  of  force  transfer  between  "soft"  and 
"hard"  material  components  such  as  layered  structures  or  fiber  reinforced 
composite  materials.  The  formulation  of  all  the.se  models  derives  from  the 
common  assumption  that  the  "soft"  material  responds  only  in  shear,  while 
the  "hard"  matter  carries  normal  loads  alone. 

When  applied  to  load  transfer  mechanisms  in  uniaxially  reinforced  poly¬ 
meric  composites  this  assumption  appears  to  provide  an  excellent  approx¬ 
imation  in  many  important  circumstances.  For  instance,  under  tension 
parallel  to  the  fiber  direction,  the  fibers  carry  between  95  and  99  percent 
of  the  applied  load.  It  can  also  be  shown  that  their  shear  distortion 
amounts  to  no  more  than  5  percent  of  their  normal  strain. 

Since  our  analysis  of  the  eccentric  case  is  constructed  as  a  perturbation 
from  the  concentric  configuration,  it  is  helpful  to  formulate  the  concentric 
case  in  a  form  which  is  readily  amenable  to  the  perturbation  scheme. 

2.1  A  Shear-Lag  Model  for  the  Concentric  Case 

Consider  the  central  fiber  of  radius  a  positioned  within  a  perfect  hexag¬ 
onal  array  of  neighboring  fibers.  Also,  let  h  denote  the  distance  between 
the  centers  of  adjacent  fibers  and  b  the  outer  radius  of  the  "shear  zone"  as 
shown  in  Figures  3a  and  3b.* 


*  Note  that  since  {albf  ^  v/  the  "shear  layer"  a  <  r  <  b  extends  bevond  the  representa¬ 
tive  volume  element  which  is  used  in  analyses  of  effective  properties.  Therefore  the 
two  cylinders  with  radii  a  and  b  should  not  be  confused  with  the  concentric  cylin¬ 
ders  employed  in  "the  three  phase"  model  (Chritensen,  1979). 
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The  fiber  volume  fraction  is  given  by  whereby 

h  =  ^  f)  CL  and  h  =  li  -  a  -  (V27i4VT  v  f)  - 1 )  u 

Let  c,  0,  z  denote  cylindrical  coordinates  and  consider  a  perfectly  hexago¬ 
nal,  uniaxially  reinforced  fibrous  array  of  infinite  extent  under  uniaxial 
stress  cTz  =  CTo.  Consider  the  case  of  a  single  broken  fiber  at  z  =  o.  The  effect 
of  this  break  can  be  analyzed  by  superposition  of  the  "undisturbed"  solu¬ 
tion  (a^  -  Go  everywhere)  and  the  solution  to  the  disturbance  caused  by 
[gz  {r,d,  o)  =  -Go  {o  <r<  a,  o<9<  Itt)  . 

The  solution  to  the  latter  problem  will  be  based  upon  the  following  as¬ 
sumptions: 

(a)  the  only  non-vanishing  displacement  is  Uz=Uz(r,z). 

(b)  discard  the  shear  deformation  in  the  broken  fiber,  whereby 


u^=z  f(z)  o<r<a  , 

0  <  Z  < 

(2) 

and 

ct4  =  E  ff'  (z) 

(3) 

As  in  the  Bernoulli-Euler  beam  theory,  the  shear  stresses  that  act 
upon  the  broken  fiber  can  be  determined  as  a  reaction,  rather  than 
from  constitutive  relations. 


discard  the  contribution  of  normal  stresses  to  the  response  of  the 
matrix  region  a  <  r  <  b.  Consequently  consider  only  the  shear  strain 

-j  m 

=  shear  stress  rz=  n  and  equilibrium  governed 


r7  r 

=  0. 
r 


The  displacements  vanishes  dX  r  =  b. 


Obviously,  it  is  necessary  to  satisfy  the  displacement  continuity  condition 
at  r  =  a  and  "global"  equilibrium  for  the  broken  fiber. 
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In  view  of  the  foregoing  assumptions  we  construct  the  solution  to  the  pre¬ 
sent  problem  as  follows; 

By  hypothesis  u(  =  fiz),  consequently  consider  uP=f(z)g{r)  [a  <  r  <  h, 

0  <  Z  <  oo). 


From  the  equilibrium  condition  (4)  we  get  =  0,  hence  g  {r)  = 

A  +  B  inr.  The  conditions  (z,  a)  =  (z,  a)  and  wf  (z,  b)  0  yield 


u 


In  {btr) 
In  (bla) 


/(z) 


(5) 


Finally,  force  equilibrium  in  the  z  direction  for  the  fiber  gives 

ica^  -  ^  ^  +  Inaxpi  (a,  z)  =  0 
dz 

which,  upon  employment  of  equations  (3)  and  (5),  gives 

F(z)-^/(z)  =  o 
where 

Ef  In  {bla) 


(6) 


(7) 


Finally,  the  boundary  conditions  cr/  =  -Oq  at  z  =  0  and  cr-^  =  0  as  z 


f(^z)  =  e'^ 

Ef  k 


and  ai=  -ctq  e'^ 

By  superposition,  the  total  stress  on  the  broken  fiber  is 


oo  yield 

(8) 

(9) 
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2 . 2  Ficcentric  Shear  Lag 


It  is  conjectured  that  a  reasonable  estimate  of  the  effects  of  random  fiber 
spacings,  exhibited  in  Figures  1  and  2,  can  be  obtained  by  analyzing  the 
idealized  circumstance  of  a  perfect  hexagonal  array  consisting  an  eccentric 
central  fiber  as  shown  in  Figure  4.  (The  circle  drawn  in  dashed  lines  shows 

the  concentric  position  of  that  fiber.)  In  the  sequel  we  let  rj  =  ejh  and  con¬ 
struct  a  shear-lag  model  which  employs  rj  as  a  perturbation  parameter. 
We  shall  assume  that  the  circumstance  of  maximum  eccentricity  within  the 
idealized  configuration  shown  in  Figure  4  (when  e  -  h  -  a,  rj  =  (b  -  a)lh)  ac¬ 
counts  for  the  effects  of  the  commonplace  randomness  of  the  arrays  shown 
in  Figures  1  and  2. 

Obviously,  when  the  inner  fiber  breaks  at  z  =  0  the  eccentricity  will  give 
rise  to  non-uniform  shear  stresses  around  the  fiber/matrix  interface. 

These  stresses  may  vary  not  only  along  the  z  direction,  but  also  with  the 
angle  6 ,  and  they  attain  their  highest  values  at  the  place  of  nearest  ap¬ 
proach.  Consequently,  that  fiber  will  be  subject  to  bending  moments  and 
undergo  rotations  in  addition  to  normal  deformation. 

In  line  with  assumptions  (a)  through  (d)  of  the  previous  section  we  con¬ 
sider  now 

u  (f,  0  ,  z  )=  f  (z  )+  T]  p  {z  )  r  cos9  (o  <r<  a,  0<P<  2k,  0  <  z  <  oo)  (10) 


Equation  (10)  implies  that  fiber  deformation  satisfies  the  hypothesis  that 
"planes  remain  planes." 


Within  the  matrix  region  ^  ^  f  ^  displacement  and  the  shear 

stress  rz  depend  on  0  and  equation  (4)  must  be  replaced  by 


m  -1 

,  1  _  0 
dr  dd 


(11) 


with 


Tr'?  =  G 


du 


m 


m 


dr 


and 


du 


m 


^0z  ~  — 


dd 


(12) 


Let  6,z  )  =  q{z)p{r,0)  then  employment  of  (11)  and  (12)  yields 
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1  ^ 

+  ' . + 


0 


vv  hereby 


I  \ 

=  p  {z)  I A  +  B  Inr  +  ^  iA^r"'  +  B,nr'"‘)  cos  md  j 

1  m=\  I 


(13 


where  even  symmetry  in  0  was  considered  without  loss  of  generality. 


Expressions  (10)  and  (13)  are  readily  amenable  to  a  solution  by  the  per¬ 
turbation  method  of  Fames  and  Beltzer  (1986),  employing  the  eccentricity 
ri  as  the  perturbation  parameter.  Although  this  method  can  be  carried  to 
any  desired  degree  of  accuracy  by  expanding  to  ever  higher  powers  of  7]  , 
we  deliberately  refrain  from  going  beyond  the  first  power  in  rj  in  deriving 
our  estimate  for  the  effects  of  fiber  eccentricity.  We  eschew  the  employ¬ 
ment  of  higher  order  expansions  because  the  idealization  involved  in  se¬ 
lecting  the  configuration  shown  in  Figure  (4)  to  represent  the  behavior  of 
random  arrays,  such  as  shown  in  Figures  1  and  2,  leads  us  to  believe  that 
higher  order  mathematical  accuracy  may,  in  fact,  be  devoid  of  physical 

content.  Nevertheless,  an  expansion  to  0(7]^)  is  given  in  Appendix  A  to 
demonstrate  the  applicability  of  the  perturbation  technique. 


Employment  of  the  expression  (10)  and  the  form  given  in  equation  (13), 
together  with  the  requirement  that  (a,  9,  z)  =  (a,  0,  z)  and  the  as¬ 

sumption  (b,  0,  z)  =  0  yield 


M  ?  (r,  9,z)= 

In  (b/a) 


In  (b/r  )  +  7]  ~  cos  6 

b^-a^ 


T]  P(z  )  — ^ -  r|  cos  6  a<r<b,  0<6<2n:  (^4^ 


b  ^  -  a  ^ 


The  normal  stress  in  the  fiber  and  the  shear  stress  in  the  matrix  are  given 
by 


/' 


(z)  +  B  b  (z)  rcos  6  J 


0<r<a 


(15) 


and 
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f"  =  m  f-J  -  ;  I  +  l|/)  .1 0),v  0 

\lnih/a>['-  2 


+  7]  /?  (:''  ~ -  1  + 

h  -  -  a  -  ^  r  ~  ' 


ros  B  ! 


!  1()  ) 


a  <  r  <  h 


l  orcc  equilibrium  for  the  fiber  in  the  z  direction  reads 


Jo  lo”^  ■‘'Jo''  ^  ^ 

dz 


!  ">  ' 
1  /  ; 


while  the  moment  equilibrium  of  the  fiber  about  the  y-axis  |0=±^- 


gives 


cos  6  drdO  + 


{2n 

io 


(a,  6,  z  )  cos  6  dB  =  0 


(18) 


upon  substitution  of  equations  (15)  and  (16)  in  expressions  (16)  and  (18) 
we  obtain  the  following  field  equations  for  f(z)  and  p(z): 


and 


"  2  -1  'yn 

f  (z)-^f(z)  =  0  ,  k^  =  — - 

Eflnibla) 


p"  {z)-J^f3(z)-^f{z)  =  0 


a 


a 


(19) 

(20) 


■?  2 

where  /  ^  +  h  n  2  _  8Gm _ t _ 

^ f  (b^  -  a^)  In  (b/a  ) 

Note  that  equation  (19)  is  identical  with  the  concentric  result  given  in 
equation  (6).  When  higher  order  expansions  in  tj  are  employed  this  is  no 
longer  the  case  and  the  field  equations  for  f(z)  and  P(z)  are  fully  coupled. 
The  uncoupling  in  equation  (19)  occurs  because  the  lowest  order  correction 

to  w  4  is  0  (ri  2) . 


For  an  extended  medium  and  a  semi-infinitely  long  broken  fiber,  0  <  z  <  oo^ 
we  have  /  (z)  =  and  the  boundary  condition  at  z  =  0 


9 


i)  lo^  if,  0,0)  nlrdO  =  -na-o  ^ 


which,  together  with  expression  (14)  gives /A  =  — ’  -  whereby 

Ef  k 

/'(::)=  ( o  j  ) 

Ef  k 

as  in  equation  (8). 

Since  the  disturbance  introduced  by  fiber  breakage  involves  no  applied 
moment  at  c  =  0,  we  have 

[()  io^  '  0,0)  cos  0  dr  dO  =  0 

whereby,  in  view  of  equation  (15),  we  obtain  the  boundary  condition 

fi'{0)  =  0  (22) 

Substituting  (21)  into  (20),  together  with  the  condition  (22)  and  the  re¬ 
quirement  that  Urn  /3  (z)  =  0  we  obtain  after  several  manipulations: 

Z  oo 


p  (z)  =  n  2  £0 

£ / 

For  realistic  values  of  V/-  (0.4  <\f<  0.7)  k  and  /  are  of  similar  magnitude. 

With  f(z)  and  /3(z)  thus  determined,  equations  (10),  (14)  -  (16),  provide  the 
"complete"  solution  to  the  perturbed  field  caused  by  the  break  of  the  inner 

fiber.  The  total  field  is  given  by  adding  cr4  =  to  expression  (15). 

We  are  now  in  a  position  to  assess  the  overload  caused  in  the  nearest 
neighboring  fiber,  centered  at  r  =  h,  0  =  0.  Assume  that  the  fiber  is  over¬ 
burdened  by  the  shear  stresses  which  are  shed  upon  it  over  an  arc  of 
length  2nbl6.  As  shown  in  Figure  5,  these  shear  stresses  correspond  to 

r  r'z  (b,  2,  &)  spanning  the  range  -  a  <  6  <  a,  where  a  =  -3-  +  a'  and 


a_ 

Ik 


(,z/a)  .  a_ 

I 


Ula) 


(23) 
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riic  overload  (force)  is  given  by 


N  (z)  =  f"  T  /?  {h,  6,  z)  hdOdz  (24) 

J  CO  J  -  ( 

Substitution  of  equations  (16),  (21)  and  (23)  into  expression  (24)  yields, 
after  several  manipulations 

N  (z)  =  a  Q  -  tj  sin  a 

+  ^-lU/a)  \ 

(a^  +  b^]c^  J( 

In  (25)  =  2  +  b^)  In  (b/a)  -[b^  -  aA 


b — +  - ^(llh -  [fi  Q-k{zla) 

(/,2.^2)c2  j 


In  addition  to  the  normal  overload,  the  shear  stresses  which  act  at  the  in¬ 
terface  between  the  matrix  and  the  nearest  neighboring  fiber  introduce 
bending  into  that  fiber.  To  assess  the  bending  effect  assume  that  these 
shear  stresses  are  distributed  uniformly  over  the  arc  PF'  shown  in  Figure 
5.  This  would  place  the  line  of  action  of  N(z)  at  the  center  of  gravity  of  the 

arc  PP',  hence  with  a  lever  arm  p  =  3VTa/2K  about  the  center  of  the  nearest 
neighboring  fiber. 


Combining  the  normal  and  bending  effects  we  have  an  overstress  on  the 
nearest  neighbor 


min  jia^  “  Ka^ 


(26) 


This  expression  indicates  that  the  overload  would  substantially  enhance 
the  tensile  stresses  at  the  point 


3. 


RF-SULTS  AND  CONCLUDING  REMARKS 


The  overloads  borne  by  the  fiber  which  is  nearest  to  the  broken  filament 
arc  shown  in  Figures  6  and  7  for  E-glass/epoxy  and  graphite/epoxy,  re¬ 
spectively.  For  E-glass/epoxy  we  took  v^=  0.65,  £/=  45  x  10^  psi  and 
G ffi  =  0.2  X  10^’  psi  ,  while  for  graphite/epoxy  we  employed  vy=  0.65, 

45  X  10^  psi  and  Gm  =  0.2  x  10^  psi. 

The  results  in  Figures  6  and  7  exhibit  the  overstresses  S  -  oJoq  at  the 
center  of  the  nearest  neighboring  fiber,  at  its  point  of  contact  with  the 
broken  filament,  P^,  as  well  as  at  the  point  which  is  diametrically  opposite 
the  point  of  contact,  P.^ .  The  total  stress  enhancements,  as  given  by  the 
factor  (S  +  1),  are  shown  by  the  curves  ac(z)/Go  and  in  figures  8 

and  9. 

Note  that  a  maximal  stress  increase  of  about  100%  occurs  in  both 
glass/epoxy  and  graphite/epoxy.  This  contrasts  with  rises  of  about  10%  - 
15  %  evaluated  by  Hedgepeth  (1961).  Furthermore,  the  effect  of  a  fiber 
break  is  confined  to  a  smaller  distance  in  glass/epoxy  than  in 
graphite/epoxy  since  the  higher  modulus  of  the  graphite  fibers  introduces 
a  stronger  stress  channelling  effect  due  to  the  magnified  anisotropy  of  the 
composite. 

An  accounting  for  effects  of  overload  on  fibers  near  the  broken  filament 
was  incorporated  by  Zweben  and  Rosen  (1970),  where  all  overloads  were 
considered  to  be  purely  tensile  and  all  fiber  arrays  were  assumed  periodic. 
They  noted  that  fibers  adjacent  to  the  broken  one  would  fail  at  the  stress 
<Tf/K,'wherQ  <jfis  the  nominal  failure  stress  and  K  denotes  the  overload 
factor.  Since  fibers  within  composites  are  known  to  break  even  under  low 
load  levels  the  consideration  of  overload  effects  appears  to  be  realistic. 
This  argument  applies  in  the  present  work  as  well,  with  the  additional 
proposition  that  stresses  due  to  fiber  bending  should  be  incorporated  in 
the  statistical  analysis.  Furthermore,  it  should  be  noted  that  due  to  the 
randomness  of  fiber  spacing  within  the  transverse  cross  section  the  factor 
K  should  be  viewed  as  a  random  variable  and  the  statistical  analysis 
should  be  reformulated  to  account  for  this  fact.  Finally,  note  that  normal 
stresses  exist  in  the  matrix  (and  fibers)  in  the  transverse  direction  in  the 
vicinity  of  a  fiber  break.  A  method  to  determine  these  stresses  was  devel¬ 
oped,  within  the  context  of  shear  lag  models,  by  Goree  and  co-workers 


(see,  e.g.,  Dharani,  L.R.,  Jones,  W.F.  and  Goree,  J.G.,  1983).  Obviously,  such 
stresses  would  cause  bending  in  neighboring  fibers  even  in  perfect  arrays. 

In  view  of  the  above  considerations,  it  is  clear  that  failure  mechanisms 
quite  distinct  from  tensile  fracture  can  occur  in  fibers  adjacent  to  a  broken 
filament.  Since  it  is  reasonable  to  assume  that  imperfections  of  several 
kinds  (such  as  weak,  kinked,  misaligned  and  broken  fibers,  debonds  along 
fiber/matrix  interfaces,  voids  and  cracks  within  the  matrix,  etc.)  are 
commonplace  within  fibrous  composites,  it  is  to  be  expected  that  the 
complex  state  of  stress  presented  in  this  work  should  be  commonly 
encountered  in  the  vicinities  of  those  imperfections.  Consequently,  the 
failure  of  filamentary  composites  may  well  be  attributed  to  complex, 
mixed-mode  mechanisms  which  are  inherently  different  form  those 
represented  by  the  strength  statistics  of  "dry"  fiber  bundles,  which  is 
associated  with  purely  tensile  response. 


Acknowledgement 

This  research  was  supported,  in  part,  by  the  Office  of  Naval  Research  under 
Contract  N00014-90-J-1556,  to  one  of  the  authors  (YW).  Dr.  Y.  Rajapakse 
of  the  Mechanics  Division,  Engineering  Sciences  Directorate,  has  been  the 
program  manager. 


APPENDIX  A:  The  Perturbation  Expansion 


Consider  two  eccentric  circles  of  radii  a  and  b,  respectively  and  with  eccen¬ 
tricity  e  as  shown  in  Figure  A-1. 

As  indicated  in  Equations  (2)  and  (5),  the  solution  to  the  concentric  shear- 
lag  problem  is  given  by 


u(  =  =/(z) 


(A-1) 


u 


m 

z 


uf{r,  1)  = 


In(blr) 

In(bla) 


/(z) 


(A-2) 


where  according  to  Equation  8 


f(z)  =  ^  ^  g-k(.z/a) 
Ef  k 


j^2  _  2Gm 

Ef  In  (b/a) 


(A-3) 


In  the  eccentric  case,  considerations  of  static  equilibrium  for  the  fiber  lead 
to  the  incorporation  of  a  bending  component  in  its  deformation.  Upon  as¬ 
suming  that  undeformed  cross-sectional  planes  of  the  fiber  remain  plane 
after  deformation,  we  choose  in  accordance  with  Equation  (10): 


ui  =  i4(r,  e,  z)=J{z)+  T]  P  (z)  rcos  0 


(A-4) 


where  rj  =  elb.  In  the  sequel  t\  serves  as  our  perturbation  parameter.* 


In  addition,  by  equation  (10)  -  (12),  equilibrium  of  the  matrix  region 
yields 


u 


m  — 
z  “ 


ln(b/r) 

ln(b/a) 


/(z)  + 


I 

-m  =  I 


COS  mO 


P(z) 


(A-5) 


It  is  now  necessary  to  satisfy  the  conditions  that  =ufz2iX  r=  a  and  that 


*  At  first  glance  it  may  appear  that  the  representation  (A-4)  is  somewhat  arbitrary. 

However,  it  can  be  shown  that  if  one  starts  with  ul{z)  =  f{z)  +  P  (z)  r  cos  6  than  a  sys- 

tematic  expansion  in  ti  indeed  yields  (z)  =  tj  /J  (z).  We  chose  (A-4)  as  our  starting  ex¬ 
pression  to  circumvent  non-essential  details. 


iit>^  =  0  on  the  boundary  of  the  eccentric  circle  of  radius  r-h. 

According  to  the  perturbation  scheme  of  Beltzer  and  Fames  (1986)  the 
value  of  a  function  F(r,  6)  on  the  boundary  Q  of  an  eccentric  circle  of  ra¬ 
dius  b  with  eccentricity  rj  =  elh  can  be  expressed  by  a  perturbation  expan¬ 
sion  in  77  as  follows 


fIc.  =f 


{o) 


-I- 


F<‘'  +  o<J., 


(o) 


ri  + 


zr(2)  ,  ^  ^ 

F  +0^1  +0*^2 


Ti^ 


+  [F®  +  o<*.r-'  +  o<J>2"'  +  o-*-3"'']r)‘  + 


(2) 


(1) 


(o) 


(A  6) 


where,  in  the  present  case 


=  -b  cos  6  Fy^ 


0^1 


b  cos^  0  fHI  -  sin^  6  f! 


and  o<^3 


^  cosd 
2 


sin^  6  F..^  -  ^  cos^  9  F, 
3 


rr 

0)  ' 

rrr 


(A-7) 


In  the  present  problem  the  function  F(o)  is  given  in  equation  (A-4),  and 
the  boundary  condition  at  u^lc„  =  o  is  satisfied  by  the  employment  of  ex¬ 
pressions  (A-6)  and  (A-7).  On  the  other  hand  the  continuity  condition 

dil  r  =  a  does  not  involve  the  perturbation  scheme  since,  by  hy¬ 
pothesis,  r  =  a  is  a  "concentric  boundary."  Obviously,  the  presence  of  two 
functions  of  z,  f(z)  and  p(z),  in  equation  (A-4)  will  necessitate  two  series 

expansions  in  (A-5)  (to  assure  r  =  a).  Denote  the  two  series  ex¬ 

pansions  by 


X  [A,nr'^  +  ^ 


Lm  =  1 

and 


m 


Lm  =  1 


COS  md 


X  [Amr'^  +  ^  C05  md 


Az) 


P(z),  respectively.^ 


Then,  a  systematic  application  of  the  conditions  «r/Co  =  0 
(a,  6,  z)  =  i4  (a,  9,  z)  yields 


(A-8) 


*  It  is  of  course  always  permissible  to  include  the  "concentric"  solution  A  +  B  Inr 
within  the  perturbation  expansion  whenever  necessary. 
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u  f  (r,  e,  z)  =  +  t;  (4 1  r  +  -  J  ^os  6 

\[n  {h/a)  '  f  ’ 

+  77^  C2ln(r/a)  +  iA2r^  +  ^\  cos  26 

[  I  /-2I  J 

+  7]^  |c3/'  r  COS  0  4- 1/4  3/'^  +  COS  36  I  /(z) 


fl2  fj^  .  f2 

b^-a^ 


cos  6 


+  772  Koln(rla)+  A2r^  +  —  cos  26 

r2  , 


+  77^  \C\r  +  ^cos6+  Azr^  +  ^  cos  36  \  Piz) 


0  (774). 


(A-9) 


where,  upon  denoting  Aq  = - - ,  A  =  — ^ — ,  we  have 

ln{bla)  b^  - 

Ai  = - 6 = 


i - —  A.  ,  B2  =  I-  -aihl 


-42=  -■! - 2 - 

M*  '  -  a^f 


(b^-a^f 


C2  = 


b^  -  log{bla) 


^3  =  -^^ - ^ - 7^0  ,  ^3  =  ^  -«-  --  Ao 

^  [b^-a'^)  3[b^-aA 


c,  = - 2 - - +  —bl - i -  Ao 

b^  -  a^\^b'^  -  a'^f  b^-  (^/«) 


D2  =  -^ - al - +_Jll - i -  Ao 
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a  II  ci 


A'. 


hA 


C^  -  - 


log  (hta) 

L  2 


h^A 


a 


,  Bi 


a- 


1  + 


log  {hta)  /,4_^4j 


D,  =-alb.L 
b~  -  a‘ 


1  + 


+ 


bl±jil 


log  {b la) 


A 


A-i  -  - 


b^  -  b"^  - 


^  +  a"*  a 


-  a*^  b^  - 


a^b'^  A 
b^  - 


The  equilibrium  of  force  and  moment  on  the  fiber,  as  given  in  equations 
(17)  and  (18),  yield  two  coupled  ordinary  differential  equations  for  f(z) 
and  )3(z).  It  turns  out  that /(z)  and  P(z)  take  now  the  following  forms 


f{z)  =  ^  +  a2T]^  +  a37]^|  +  (v277^  +  V377^)  e'^(2/a) 

\Efk  I 

(A-10) 

^  A 

Piz)  =  [Piri  +  P2'n^  +  )3377^)  +  (yjTy  +  72??^  +  73^7^) 


/s 

with  k  =  k  +  r\k\  +  ri'^ki  +  rj^k^  /  =  /  +  77/1  +  77^/2  +  77^/3. 


As  noted  earlier,  we  refrain  from  employing  the  higher  order  results  of 
this  Appendix  in  our  computations,  because  the  approximations  involved 
in  our  model  do  not  warrant  their  use.  In  fact,  the  utilization  of  higher  ex¬ 
pansions  may  convey  the  fallacious  impression  of  higher  accuracy. 
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Figure  Titles: 


Fig.  la:  A  typical  cross-section  of  graphite/BMI  composite,  estimated  vol¬ 
ume  fraction  vf  =  55%  (courtesy  of  Professor  W.W.  Stinchcomb)  (xlOOO). 

Fig  lb:  A  typical  cross-section  of  graphite/epoxy  composite,  estimated 
volume  fraction  vf  =  65%  (courtesy  of  Professor  W.W.  Stinchcomb) 
(XlOOO). 

Fig.  2a:  Monte-Carlo  simulation  of  a  cross-section  with  fiber  volume  frac¬ 
tion  Vf  =  40%. 

Fig.  2b:  Monte-Carlo  simulation  of  a  cross-section  with  fiber  volume  frac¬ 
tion  Vf  =  65%. 


Fig.  3:  A  perfect  hexagonal  fibrous  array,  with  broken  central  fiber^  sur¬ 
rounded  by  a  concentric  "sleeve"  of  matrix  material  in  shear. 

(a)  Cross  sectional  view,  (b)  Side  view. 

Fig.  4:  A  hexagonal  fibrous  array  as  in  figure  3,  with  broken  eccentric  fiber 
eccentricity  e,  and  nearest  neighboring  fiber 

Fig.  5:  Geometric  considerations  for  evaluating  the  overload  carried  by  the 
nearest  neighboring  fiber  @  due  to  the  fracture  of  the  fibei^^.  The  load 
shed  upon  is  related  to  the  angle  a  (see  text). 

Fig.  6:  The  dimensionless  overload  stresses  on  the  fiber  S(Pcz),  S(Pq,z) 
and  S(P/{,z)  at  points  P^,  Pq,  and  P/\  vs.  the  non-dimensional  distance 
z/a.  Case  of  glass/epoxy,  with  vf  =  50%. 

Fig.  7:  Same  as  figure  6,  but  for  the  case  of  graphite/epoxy  and  vf  =  65%. 

Fig.  8:  The  dimensionless  stresses  Oc{z)/g^  and  (t^(z)/<To  carried  by  the 
fiber@at  points  Pq  and  P\  vs.  the  non-dimensional  distance  zta.  Case 
of  glass/epoxy,  with  vf  =  50%. 

Fig.  9:  Same  as  figure  8,  but  for  the  case  of  graphite/epoxy  and  vf  =  65%. 
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